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ABSTRACT 

We present the full tree-order calculation of the process e + e~ — > li>i qq' 
including the width effects in a way which is consistent with gauge in- 
variance and high-energy unitarity and which has a natural interpretation 
within standard perturbation theory. We also include C and P preserv- 
ing deviations from Standard Model trilinear vector boson couplings. An 
analysis of the sensitivity of these channels on the trilinear couplings, us- 
ing well measured angular variables of the four-fermion final state, is also 
presented. 



February 1995 



LEP II will provide us, for the first time, with the opportunity to study directly 
the W boson, above the energy threshold for IF-pair production. One of the most 
important measurements at these energies will be the determination of the trilinear 
vector boson couplings [I, 2, 3, 4, 5], a characteristic manifestation of the underlying 
non-Abelian symmetry of elementary particle interactions [6]. In order to determine 
these couplings we need accurate predictions of the cross sections which will be mea- 
sured at LEP II. It is therefore evident that a good theoretical description of the four 
fermion production is indispensable. 

In this paper we present a full tree-order calculation of the processes: 



taking into account the non-zero width of the heavy unstable particles W and Z 
in a way which is consistent with gauge invariance and well justified within standard 
perturbation theory. We also include all possible deviations 1 from the Standard Model 
trilinear vector boson interactions which respect C and P invariance. The need for this 
calculation stems from two different directions: from the experimental point of view it 
is evident that the LEP II environment will consist of four-fermion final states; from 
the theoretical point of view, we expect that finite width effects as well as non-resonant 
contributions, estimated to be of the order of a few per cent, will be significant for 
the trilinear vector boson couplings measurements. 

The Feynman diagrams contributing to the processes Eq.(I) are given in Figs. 1, 2. 
Graphs I to 3 contribute to IF-pair production, whereas graphs I to 6 as well as II to 
16 contribute to single- IF production, e + e~ — > e~v e W + [4, 5]. Graphs I, 2, II and 12 
are the signal graphs, since they include the trilinear boson vertex function, whereas 
graphs 1 to 3 are the double-resonant graphs, 4 to 16 are the single-resonant ones and 
finally graphs 17 to 20 are the non-resonant ones. 

The amplitudes have been calculated using the E'-vector formulation [4] of the 
spinor product technique [7]. This is based on the following representation for the 
fermion 'current', for massless fermions: 
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Sometimes referred to as anomalous couplings. 
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+ + (Plx + iPly)(P2x - iP2y) 



E- = VP1P2 ~ K1 ^ (3) 



, pf pi 

with p^ = p° ± p 3 . 

The main advantage of this formulation is that there is a very straightforward way 
to write down the expression for a given Feynman graph. For example, for the graph 
number 7 the expression is simply (omitting propagators and coupling constants): 

T~ Yl sign(pi)E-(p q ,pi) ■ Ex(p e +,p e -) E-(pi,pq,) ■ E-(p t ,p Vt ) (4) 

i=q r e~ ,e+ 

where sign(p) = ±1 for outgoing and incoming momenta respectively. Furthermore 
the FORTRAN code representation of the amplitudes is very simple, facilitating consid- 
erably the writing and checking of the programme. 

One of the first problems, inherent in any four-fermion calculation, is the treatment 
of the propagator of the heavy unstable particles, W and Z, appearing as intermediate 
states in all Feynman graphs [8]. The problem is that if we implement the tree-order 
form of this propagator 
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we have a non-integrable singularity in phase space, located at p 2 = m\. The well- 
established procedure to cure this problem is to take into account the finite width of 
these particles. The propagator takes the well-known Breit-Wigner form 

f 



p 2 — m v + imyVvyp ) 

From the point of view of the quantum held theory, the Breit-Wigner propagator can 
be understood as the effect of the resummation of all self-energy graphs as shown in 
Fig. 3. The 'width' of the particle is then identified with the imaginary part of the 
one-loop self-energy contribution and it is a function of p 2 . The 'fixed' width Ty 
is now defined at p 2 = rriy, and it is connected to the decay of the on-shell heavy 
particle, as is dictated by the optical theorem. 

Although all these seem to be rather elementary, the treatment of the width effects, 
using the above prescription, leads in general to gauge-invariance violation [9], as well 
as to unitarity violation at high energies. Let us hrst see how this problem appears 
in the 't-channel' graphs ff, 14, 17 and 19 of Fig. 2, where a photon is exchanged 
carrying a space-like momentum. Their expressions, omitting trivial factors, are given 
as follows: 



1 1 1 



(P1-P3) 2 (P4 + P5e) 2 2p 5 -p 6 -m 2 v 

u 3 7^i ^57^6 ^27^(^4 + ^56)7^4 
1 1 1 



14 



(.P1-P3) 2 -2p 2 -p4-m 2 v 2p 5 -p 6 -m 2 v 



2 



lt 3 7^1 U 5 -f a U 6 U 2 Y U 4 



1 1 1 



7 



19 



(P1-P3) 2 (P1-P3-P5) 2 -2p 2 -p 4 -m 2 y 

M 3 7^Ml U 5 7^Gfo + ^5 - Pl)7^6 U 2 ^U A 
1 1 1 

(P1-P3) 2 (pi-P3-Pe) 2 -2p 2 -p 4 -m 2 v 

U 3 7^1 M 5 7^(i>l - i>3 - k,)lvU 6 U 2 Y U 4 (5) 

where Pi, P3, P4, Ps, P6 are the momenta of the incoming electron and the outgoing 
electron, neutrino, quark and antiquark, respectively, and T 1 is the three-boson 
vertex function. The singularity appearing in the collinear limit Ps||pi is integrable due 
to the non-zero electron mass and the result is nothing but the well-known logarithmic 
rise of the total cross section, log(s/m 2 ), where m e is the electron mass, which forms 
the basis of the Weiszaker- Williams approximation [10]. If we look more closely at 
the expressions given above, we see that in order to get a logarithmic behaviour (or a 
logarithmic divergence in the case of massless electrons) we need to have the following 
relation 

(px - p 3 ) 2 (T n + T 14 + T 17 + T 19 ) = 

in the limit P3 = x p^, which is nothing but the gauge invariance at the subprocess 
level 7 + e + — > v e qq' . Let us now write the above equation in the collinear limit: 
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(P4 + P5 + Pe) 2 2p 2 • p 4 + m 2 w -18) (p 5 + p 6 ) 2 - ra 2 w + 
ST^E^E^qPi 1 E 56 ■ E 2 a 



(ps + Pe) 2 - rn'w + 27 1 - x 2p 2 • p 4 + m 2 ^ - iS 

where we have allowed for several possible contributions, including the width of the W 
both in time-like and space-like regions (7 and 8 respectively) as well as a modification 
of the three-boson vertex. The solutions are easily worked out and they can be 
enumerated as follows (ST /lafjE^E^pi = 8T E^q ■ £"24/(1 — x)): 

1. <5 = 7 = (T = which is the standard tree-order result. This of course has the 
disadvantage that the matrix element squared is not integrable over the phase 
space. A 'solution' to this problem, hrst proposed in reference [11], is to multiply 
the tree-order amplitude by a factor 

2 2 
q - m w 



q 2 — m w + imw^w 

where q 2 is the (positive) momentum squared of the intermediate W. This of 
course provides a gauge-invariant quantity, as an approximation of the tree- 
order cross section, which, nevertheless, is not well justified within standard 
perturbation theory. 
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2. 6 = 7 = rriw^w and <5r = which is the so-called 'naive scheme', where a 
fixed width is considered in any W propagator. This 'solution' shares the same 
problem as the first one, since there is no justification of the width-term in 
propagator functions with space-like momentum flow, where the imaginary part 
of the self-energy contribution is vanishing. 

3. Finally, the last solution is given by: 

7 = ST = ^-^(p 5 + Pef- 
mw 

It is exactly this last solution which is the basis of our approach to the width problem. 
What it suggests is that we have to take into account the relevant contribution from 
the one-loop correction to the three-boson vertex [12]. This is shown in Fig. 4, and 
we will use the term Breit-Wigner triangle, in order to emphasize the relation of this 
correction to the resummed Breit-Wigner propagator. 

It is very easy to verify that the imaginary part of the one-loop graph, shown in 
Fig. 5, fits perfectly the requirement of our third solution. It is not difficult to realize 
that the relation derived in the collinear limit is nothing but the standard Ward 
identity, which relates the imaginary part of the one-loop vertex with the inverse of 
the Breit-Wigner propagator 2 . 

Let us now come to the high-energy unitarity, which is related to the double- 
resonant graphs 1, 2 and 3. In the high energy limit the dominant part of each 
amplitude can be obtained by the replacements: 

— a a i a i — B [3.0 

M57 u & — > P 5 + P 6 an d u 3 -f' u 4 — >■ p' 3 + p' 4 . 
Then, all amplitudes simplify a lot and become proportional to 

■S E 2 1 ■ (P5 + P6 - P3 - Pi) , 

where s = (p\ + P2) 2 . The cancellation of these terms, whose contribution to the cross 
section grows linearly with s, is guaranteed by the following equations (g^ is the 
coupling of Z with the fermion / and gzww is the coupling with the W boson): 

Qe + g^gzww = and Q e + g^T- ] gzww + —. j = 

a consequence of the ,9/7(2) X £7(1) gauge symmetry. It is more convenient to rewrite 
them in the following form 

- - c 1 

e\ ■ 9w + = 

2 For a comprehensive analysis on the treatment of the width of unstable particles in perturbation 
theory, see reference [13]. 
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where 

eA = { ^ ) 

and 

9w = ( J 

\ gzww j 

When we consider the width effects we have to resum the self-energy contributions 
for the neutral boson propagator. The building block of this resummation is, as usual, 
the imaginary part of the one-loop contribution, which can be written as 

r-- ie *J_ N ( E f Qf(gl + g f R ) \ 

24tt 9 \ Y.fQf(gl + g f R) E f ((gl) 2 + (g f R ) 2 ) ) 

where N g is the number of fermion generations, and the summation over the fermions, 
J2f, is restricted within a complete generation. In deriving this formula we have 
neglected all fermion masses, which are irrelevant in the high energy limit. It is a 
matter of algebra to realize that the vector g*w is an eigenvector of the matrix K, and 
therefore the resummed result can be written in the form 

(Q e + gi +) gzww)(l + !e 2 ^- Ng^-)- 1 

and 12 1 

(Q e + g^gzwwKl + «e 2 - N g — y 1 + — 

Z47T S w ZS W 

for right- and left-handed electrons respectively. As is easily seen, the right-handed 
amplitude automatically satisfies the high-energy unitarity condition, since it is pro- 
portional to the tree order result. For the left-handed part we have to add the relevant 
Breit-Wigner triangle contribution. There are three cut graphs contributing in this 
kinematical configuration, but only the cut on the neutral boson site is relevant for 
the high-energy limit. In this limit the effect of the triangle graph is given by the 
tree-order result multiplied by the factor (I + i e 2 ^- iV„-§-), which cancels exactly 
the self-energy resummed contribution. This result also relies on the gauge symmetry 
relations: 

J2'Qf = ^N g and Y.'9 ( f ] = ^ g g zww 
J J 

where J2f ' means that t = 1/2 and t = —1/2 fermions are taken with a relative minus 
sign, which gives rise to a complete factorization of the Breit-Wigner triangle graph 
contribution, in the high energy limit. 

Finally we would like to briefly comment on the reason why the sum of the Breit- 
Wigner propagator and the Breit-Wigner triangle form a subset of gauge-invariant 
graphs. As far as the high-energy limit is concerned, the result is proportional to N g a, 
where a = e 2 /47r is the electromagnetic coupling constant, since the loop summation 
is taken over complete fermion generations. It is rather elementary to realize that 
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in the limit N g — > oo and a — > 0, with A 3 a=fixed, the only graphs that survive are 
exactly those considered in our discussion and therefore they are by themselves gauge 
invariant. 

It is worth while to notice that, at the amplitude level, the cancellation of the most 
divergent term, which is proportional to s/m^, holds for each fermion separately. The 
reason is that the triangle contribution is proportional to the vector 



Therefore the triangle contribution cancels exactly the resummed self-energy one, and 
terms growing linearly with s are absent, as unitarity requires. 

The bosonic contribution to the width effects should be treated more carefully. 
First of all the bosonic self-energy contribution can be resummed in the same way 
as the fermionic one and it is easily proved that the vector gw is still an eigenvector 
of the corresponding /C-matrix, in the high energy limit. Nevertheless, including the 
triangle graph contribution does not, generically, lead to a gauge invariant result 
and one has to take into account box graphs as well [13, 14]. Furthermore it is not 
possible to isolate a part of the higher order corrections, which is gauge invariant, in 
a way similar to the one suggested by the high A^-limit in the fermionic sector. A 
possible alternative is to use background held techniques [15] in order to restore the 
gauge invariance. It should be emphasized, however, that, due to threshold effects 
arising from the phase-space integration, the bosonic graphs do not contribute to the 
definition of the fixed width of the W and Z, at this order of perturbation theory. The 
same is still true for the U{1) gauge-invariance restoration in the small momentum 
transferred limit, — > 0. Finally their contribution at LEP II energies is supressed 
by factors of the form I — Arriy/s, where ray refers to the mass of W, Z and Higgs 
particle. 

After completing the amplitude calculation we now proceed with the integration 
over the phase space. To this end we have implemented a Monte Carlo integration 
algorithm which is essentially identical to the multichannel approach of references 
[16, 17]. The problem is that the amplitude we have to integrate over is a very 
complicated function of the kinematical variables, peaking at different regions of phase 
space. The idea is to define different kinematical mappings, corresponding to different 
peaking structures of the amplitude and then use an optimization procedure to adjust 




whereas the corresponding /C-matrix is proportional to 




with the property that 



Mg w = -g'w 
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the percentage of the generated phase-space points, according to any specific mapping, 
in such a way that the total error is minimized. Notice that the four particle phase 
space is determined by eight independent variables. For each algorithm, the set of 
these variables and the corresponding mapping to the eight pseudorandom numbers, 
Pi, i = 1 , . . . , 8, is given below: 

m\tf = m 2 w + m w ^w tan(p 1 (y + - y_) + y_) 

_ 1 E 2 -m 2 w 
y + = tan (- 



tan 1 ( 



rnw^w 
raw 



Tw 

m 2 l v i = m w + m w Fwt&n(p 2 (y' + - y~) + y~) 

y = tam -, (^-^-^ ) 

mw± w 

COS 6w = 2/J 3 — f , (f>w = 27T/J 4 

cos 0* = 2/?5 — 1, cf)* = 2tt pq (in the £, vt rest frame) 

cos 0* = 2pr — 1 , </>* = 27T/? 8 (in the q,q' rest frame) (7) 

which simulates the peaking structure of graphs f and 2, 

m 2 -, = + m^r^tan( / oi(y + - y_) + y_) 



9,9 

(p 2 ^ + (l-P2)m^) J 

1/(1-,) 



COS ^ = 1 - ^(1 - p 3 )(l - C mm ) (1 ^ + p 3 (l - CmaxY 1 

(f)i = 2irp 4 

cos 6p = 2ps — f, <f>p = 2tt pq (in the i>i,q,q' rest frame) 
cos 0* = 2p7 — f, </>* = 2tt ps (in the rest frame) (8) 

which for x = 1 simulates the peaking structure of graph ff and with x — > the 
peaking structure of graphs f 4, and finally 

2 2 

m - 



m w + m w T w t&n(p 1 (y + - y_) + y_) 
m\ Vt = m 2 w + m w ^w t&n(p 2 (y' + - y~) + y~) 



9,9 

2 



COS t>w 



a 



Ow = a-^l-p^a-lY^ + psia + lY 1 -^ 



i/(i-,) 



Ep w ~ ™\,v t 



E\p w \ 
4>w = 2irp4 

cos 0\ = 2p 5 — 1 , <j>£ = 2tt p & (in the I, vi rest frame) 

cos 0* = 2p7 — 1, (f>* = 2tt ps (in the q, q' rest frame) (9) 

which simulates the peaking structure of graph 3. 
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In order to study the sensitivity of the processes under consideration on the trilin- 
ear boson couplings we need a parametrization of these interactions that goes beyond 
the Standard Model . There are of course several possibilities, but we are going to 
restrict ourselves to the most economical one. The relevant interaction Lagrangian 
can be written in the following form [2] : 



£ = \ 9v(V^W~ tlv W +v - V^W+^W-" + kvV^W^W 



+ ^lk v »' w+f " w; * 



(10) 



where 



W, 



d v W; 



W 



is the VK-boson held, and g 1 = e, gz = ectg 6 W , n 1 = nz = 1 and A 7 = Xz = 
at tree order in the Standard Model . It is more convenient to express the different 
couplings in terms of their deviations from the Standard Model values. For this we 
define the following deviation parameters: 



6; 



gz 



ctg 6 h 

I 



x z = (k z - l)(ctg 6 W + 8 Z ) 



It is worth while to note that the interaction Lagrangian becomes linear with respect 
to the above parameters (including also A 7 and Xz)- 

The VW + W~ vertex function can be derived from the above interaction La- 
grangian and has the following form [4]: 



(k + X— r )q - p + 
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+ 
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9Pn 
X 
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p\ 



(k + X—^-)q 



m 2 w 



X 



+ [gap 



x q -^)( P+ 



m 2 w 



2~q ga[3P-f, 



w 



m 2 w 



X 



m 2 w 



P-p(g^aP- • q ~ P-»q a ) 



X 



m w 



q qy.g a (i 



(12) 



where p± are the momenta of respectively, and q the momentum of the neutral 
gauge boson (7, Z). Note that in many processes a great simplification occurs due 
to fermion current conservation; for example the term proportional to q^ vanishes 
whenever the neutral boson is coupled to a massless fermion pair. 

During the last few years, considerable progress has been achieved concerning the 
understanding of the physics underlying the non-standard boson self-couplings. As 
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Gounaris and Renard [18] showed, the deviations from the Standard Model couplings 
can be parametrized in a manifestly gauge-invariant (but still non-renormalizable) 
way, by considering gauge-invariant operators involving higher-dimensional interac- 
tions among gauge bosons and Higgs held. Restricting ourselves to SU(2)l X Z7(l)y- 
invariant operators with dimension up to 8, we can have the following list [5, 19]: 

Ow* = (D^Y t-W^(D^) 

Ow = ^(W^ x W" v ) ■ 
0' w „ = (&t-W^<$>)(D^)\DA) 

o' w = (r p xr,)-($tI$)B; (13) 

where t; = |<7 4 - (<t 4 - are the Pauli matrices), 

= d il B v — djjB^ 

where is the U(1)y gauge held, 

= d,w v - d v w, - g 2 W, x w v 

where W are the SU(2)l gauge helds and 

$ = ( + # + **°) ) 

is the Higgs doublet. The covariant derivative D ^ is given, as usual, by 

Dp = + i g 2 r ■ Wp - % g x Bp 

and e = g 2 sin 6 W = g\ cos 6 W . 

The interaction Lagrangian can be written now as 

C = C SM + Y.3 l O l (14) 

where Oi are the operators given in Eq.(13). The effect of each operator on the 
trilinear coupling parametrization, given in Eq.(lO), is described in table 1. The main 
consequence of the above considerations is that SU(2)l x U(1)y gauge invariance 
implies relations between deviation parameters, which must be taken into account in 
a phenomenological analysis. Furthermore, taking into account the limitting statistics 
expected at LEP II this is a very welcome result, since an investigation of the full 
multiparameter space seems to be rather impossible. In this paper, adopting this 
point of view, we focus our analysis on the six one-parameter cases, as given in table 
1. 
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case 


operator 


non- vanishing deviations 


1 




tZ/^/ • Xj 7/ tZ/^/^J ^ 1 n 


2 


Uw* — Ob* 




3 


Ow§ 




4 




A 7 , = A 7 


5 


D' 


^ 7 i X ^ ^'yCt^ $ yj 


6 


O'w 


^7, = — A 7 tg 2 ^u, 



Table 1: Contribution of different gauge-invariant operators to the anomalous cou- 
plings (s w = sin^ and c w = cos6 w ). 



Having clarified the amplitude as well as the phase-space calculation we now pro- 
ceed with the presentation of our numerical results. Several checks have been made in 
order to guarantee their correctness. First of all both the £7(1) gauge invariance and 
the unitary high-energy behaviour of the calculation have been successfully tested at 
the numerical level. Furthermore we have checked our programme with the EXCALIBUR 
Monte Carlo [16], and taking identical conditions (width implementation, input pa- 
rameters, kinematical cuts, etc) we found complete agreement at all energies at which 
we run the programme (175, 190, 200 and 500 GeV). It should be mentioned however 
that in what follows we have implemented our novel approach to the width effects, 
which is different from that of EXCALIBUR. Finally, since this calculation is an expan- 
sion of the single- W production programme [4], several checks already done for this 
process are still in order for the four fermion calculation. 

The input Standard Model parameters we have used follows: 

Mw = 80.5 GeV, Mz = 91.19 GeV, sin 2 W = 0.23, a = 1/129 

whereas the phase-space cuts are defined as: 

cos0* < 0.95 cosset < 0.95 cos 8if < 0.95 
E t ,E jet > 5 GeV and m q ^ > 10 GeV (15) 

The 'width' of the W is calculated by the formula: 

Tw ct 
mw 12s^ 

where rif = 9(12) for y/s < m top + nib(\/s > mt op + ra&). The 7 — Z propagator is 
defined by 

v = -lyi-zcDo)" 1 

s 

where 
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As we have already pointed out, lor the analysis we are attempting in this paper we 
will restrict ourselves to the one-parameter cases described in table 1. This simplifies 
the analysis, since the results (total cross sections and angular distributions) can 
always be written in the following simple form: 

(Ti = of + xa\ + x 2 a\ (16) 

where x stands for any deviation parameter (Sz, x 1 or A 7 ), cr° is the Standard Model 
result, o\ is the interference term and a\ is the signal term squared. The subscript 
i refers to a given phase-space configuration, or to a given phase-space variable. We 
find it more convenient to work with this parametrization rather than to specify some 
particular (arbitrary) values for the deviation parameters. In any case, the values of 
cr°, a\ and a\ completely determine the dependence on x. 

In Fig. 6 we show the total cross section as a function of the energy from 175 GeV 
up to 500 GeV. The Standard Model cross section reaches its maximum around 200 
GeV, whereas the interference term is negative and small in absolute value in the full 
energy range up to 500 GeV. The main contribution, especially at higher energies 
comes from <r s , and it is more than an order of magnitude larger that the Standard 
Model prediction at 500 GeV. It is worth while to mention the difference in the energy 
growth between cases 3 and 4: this is to be expected, since the relations among the 
different deviation parameters given in case 3 are expected to result in a cancellation 
of the most divergent part of the amplitude[2]. 

In order to further analyse the sensitivity of the cross section to the trilinear 
couplings, we have to use differential distributions. There are in general two different 
and in some way complementary strategies in this analysis. The 'optimistic' one is 
to use the fully differential cross section, a function of the eight independent phase- 
space variables [2, 20]. This is of course the 'ultimate' information one can get, 
but since experimentally some of the phase-pace variables are poorly measured, one 
has to cope with a low efficiency factor. The 'realistic' one is to focus on these 
kinematical variables and the corresponding differential distributions, which are well 
measured experimentally and try to maximize the sensitivity with respect to the 
trilinear couplings. In this paper we are following the latter approach and we calculate 
the differential distribution with respect to the angle of the produced lepton in the 
laboratory frame. In order to compare the sensitivity of the lepton angle with the 
sensitivity of other variables, which are a priori more sensitive to the trilinear boson 
couplings, we calculate also the following differential distributions: 

1. da /d cos 6w, where 6w is defined as the angle of the vector p q qi with the electron 
beam in the laboratory frame, which for the signal graphs is identical to the angle 
of W + , and 

2. da /d cos 6j, where the angle 6* is defined as the angle between the spatial mo- 
mentum of the outgoing quark in the rest frame of the q,q' system, and the 
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momentum p q qi in the laboratory frame, which for on-shell W decay has been 
used polarization analyser [2]. 

The Standard Model distributions are given in Fig. 7 for y/s = 175 GeV and 205 
GeV. As is evident from this figure, the distribution of the electron is very similar to 
the distribution of the W~ (which is taken from the W + distribution by the replace- 
ment cos 6 — >■ — cos 6) and peaks in the direction parallel to the electron b earn, as a 
combined effect of the neutrino as well as photon t-channel exchange contributions. 

The effect of the deviation terms on these distributions can be traced by studying 
the ratios 

^ aIld ^ 

as functions of the corresponding angular variables. As an example we present, in 
Fig. 8, the above-mentioned quantities at y/s = 205 GeV for the electron channel. 
Similar results are obtained at y/s = 175 GeV as well as for the 'muon' channel, e.g. 
when a lepton different from the electron is produced in the final state. The main 
conclusion is that the angular dependence of the deviation terms is rather important 
and the general tendency is to flatten out the background distributions in all cases. 



Operator 


0t 




9w 


X>y 


-.37 .83 


-.62 1.13 


-.46 1.06 


-.17 .35 


-.30 .81 


-.21 .86 


Ow§ — Ob§ 

Sz 


-.67 .62 


-.70 .70 


-.22 .24 


-.33 .41 


-.37 .55 


-.17 .26 


X>y 


-.24 .31 


-.26 .33 


-.08 .09 


-.11 .28 


-.13 .28 


-.06 .07 


Ow 

A-y 


-.38 .40 


-.45 .46 


-.14 .16 


-.16 .27 


-.21 .29 


-.09 .13 


0' W q> 

X>y 


-.31 .41 


-.39 .60 


-.14 .15 


-.13 .22 


-.18 .50 


-.09 .11 


O'w 


-.44 .91 


-.62 .94 


-.42 .76 


-.20 .53 


-.31 .53 


-.21 .50 



Table 2: Limits on deviation parameters from electron channel. The hrst row corre- 
sponds to y/s = 175GeV and the second one to y/s = 205GeV. 



In order to determine the sensitivity to the trilinear vector boson couplings we 
implement a standard, though idealized, \ 2 analysis. More precisely we define a '% 2 ' 
function as follows: 



x 2 = E 



si 
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Operator 




0* 

j 


0\y 


i~J R(T) 
D Sf 

X>y 


-.29 1.49 


— .45 1.60 


-.28 .38 


-.13 .19 


-.20 .95 


-.13 .17 




-.51 .56 


— .56 .65 


-.15 .16 


-.26 .38 


-.31 .46 


-.11 .14 


LJ T/T/(T) 


-.17 .30 


-.20 .35 


-.05 .06 


-.08 .17 


-.10 .28 


-.04 .04 




-.24 .36 


-.33 .46 


-.09 .10 


-.12 .21 


-.16 .26 


-.06 .07 


tC^y 


-.23 .35 


-.29 .72 


-.09 .09 


-.09 .13 


-.13 .51 


-.05 .06 


A-y 


-.31 .98 


-.48 1.00 


-.30 .50 


-.15 .51 


-.23 .52 


-.15 .32 



Table 3: Limits on deviation parameters from 'muon' channel. The hrst row corre- 
sponds to y/s = 175GeV and the second one to a/s = 205GeV. 

where e 4 - is the statistical error of the z-th bin. In the present analysis we have chosen 
six equidistant bins in the range —0.95 < cos 6 < 0.95, where 6 stands for any of 
the angles under consideration. Furthermore, the integrated luminosity at \Js = 175 
GeV is taken to be L = 500pb _1 , whereas at y/s = 205 GeV we have used the value 
L = 300pb . The errors on the deviation parameters are now given by standard 
statistical theory [21] and the two standard deviation results are presented in tables ,. 
It should be mentioned, however, that the errors given in the above tables are in some 
sense a lower limit on the corresponding 'real' ones, since they do not take into account 
the goodness of the fit statistics, which relies on how well the real data are described by 
the Standard Model predictions. Nevertheless assuming that the real data would be 
described by the Standard Model and that x 2 /dof ~ 1 we have a rather good estimate 
of the expected sensitivity. The main point underlined by this analysis is that, hrst 
of all, in most cases, the 6w is more sensitive than the other variables, whereas 6 e 
is more efficient than 6*. Secondly, the sensitivity increases with increasing energy, 
even though the luminosity is substantially decreasing. Finally, electron and 'muon' 
channels seem to give comparable results, even though in the 'muon' channel we have 
doubled the corresponding statistics by including both // and r events, which share 
the same theoretical description. The limits obtained on the deviation parameters 
are of the same order of magnitude as those extracted from on-shell W production, 
using the same angular variables. It is expected, however, that a more sophisticated 
analysis, using the method proposed in reference [20] could substantially increase the 
sensitivity to the trilinear boson couplings. Such an analysis seems to be suitable for 
the four-fermion production and it is important to be carried out, taking into account 
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the precise experimental conditions. 

In summary, we have presented a full tree-order calculation of the four-fermion 
production processes 

e + e~ — > £u£ qq . 

We have shown that a consistent treatment of the width effects is possible, at least for 
the fermionic contribution, and we have implemented it in a Monte Carlo programme. 
The C and P preserving deviations from the trilinear Standard Model couplings have 
been taken into account, and a hrst analysis on the sensitivity of the four-fermion 
production on these deviation parameters is presented. 
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Figure 2: Feynman graphs contributing to e + e — > ev e qq' in the t-channel. 
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Figure 3: The Breit-Wigner propagator 




Figure 5: The cut graph contributing to the restoration of £7(1) gauge invariance. 
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Figure 6: The total cross section, in picobarns, as a function of the cent re- of- mass 
energy for cases 3 and 4 as explained in the text. The solid, dashed and dotted lines 
correspond to the Standard Model interference and signal terms respectively. 
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Figure 7: The Standard Model distributions da /d cos 6 e (solid line), da /d cos 6w (dot- 
ted line) and da/dcos6* (dashed line), in picobarns. 
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